Parallel and series nucleation 



A note on the nucleation with multiple steps: 
Parallel and series nucleation 
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Parallel and series nucleation are the basic elements of the complex nucleation process when two saddle 
points exist on the free-energy landscape. It is pointed out that the nucleation rates follow formulas similar 
to those of parallel and series connection of resistors or conductors in an electric circuit. Necessary formulas 
to calculate individual nucleation rates at the saddle points and the total nucleation rate are summarized and 
the extension to the more complex nucleation process is suggested. 
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I. INTRODUCTION 

Nucleation is a very basic phenomenon which plays a 
vital role in various materials processing ranging from 
steel production to food beverage industries 1 -. Recently, 
researchers' interests have turned to complex material 
nucleations which are relevant to our daily life 1 -. The 
nucleation of such complex materials can also be com- 
plex and may consist of multiple steps 2 -! 3 -. These multiple 
steps correspond to the multiple saddle points on the free 
energy landscape and can occur in series or in parallel. 
Therefore, the nucleation with multiple steps can be akin 
to the text-book problem of electric circuitry. 

The most basic unit of electric circuitry is the paral- 
lel or series connections of two resistors. Likewise, there 
are some problems of nucleation where there are two sad- 
dles points on the free-energy landscape. For example, 
Ray et ali^ found the two saddle points (double bar- 
riers) on the two-dimensional free-energy landscape of 
a droplet condensation of binary vapor-. A sequential 
two saddle points on the free-energy landscape along a 
fictitious one-dimensional reaction coordinate has been 
studied by Valencia and LipowskyS. to study the con- 
densation on a patterned substrate. These free-energy 
landscapes imply two successive series nucleations. A 
similar double-barrier model was assumed to study the 
nucleation and growth by various authors 2 -^^ when an 
intermediate metastable phase exists. 

Recently, the present author— re-examined the free- 
energy landscape using capillarity approximation when 
an intermediate metastable phase existed. The location 
of the two saddle points on the free-energy landscape 
seems to suggest a parallel nucleation in contrast to the 
assumption of a series nucleation of previous authors 2 -!^. 
Figure [1] shows the free-energy landscape when two sad- 
dle points have the same free-energy barrier 8 . The two 
red points indicate two saddle points and the red line in- 
dicates the series nucleation route assumed by previous 



author a 2 ' 6 ' . This figure clearly suggests a parallel nucle- 
ation with two nucleation routes which pass through two 
different saddle points. Similar free-energy landscapes 
with two types of critical nucleus of different composi- 
tions have been predicted by Ray et ali^ for partially mis- 
cible binary systems such as diethylene glycol/ n- heptane, 
methanol/cyclohexane and dicthlcnc glycol/benzen from 
the free-energy calculation using the capillarity approxi- 
mation, and by Chen et aliS for water/ n-nonane mixtures 
using Monte Carlo algorithm with the umbrella sam- 
pling technique. Experimental result for the water/n- 
nonane mixture by Wagner and Strej*i£ is also consis- 
tent to the picture of parallel nucleation with two saddle 
points. Therefore, parallel and series nucleation seem 
to be common basic units of various complex nucleation 
phenomena. 

In this report, we consider the most basic elements of 
complex nucleation. A series and a parallel nucleation, 
which are akin to the basic units of electric circuit: series 
and parallel connection of resistors or conductors. We 
summarize the basic ingredients of the parallel and the 
series nucleations and emphasize the complex nucleation 
process which could be sorted out by mapping the com- 
plex nucleation process into a simple electric circuitry. 



II. CONTINUUM DESCRIPTION OF NUCLEATION 
PROCESS 

Nucleation kinetics can be studied using the discrete 
or continuum description o 1 ' 11 . In this report, we use the 
continuum description since we want to make correspon- 
dence to the classical electrostatics more transparent. 
In the continuum description, the discrete Master equa- 
tion is transformed into the Fokker-Planck equatio n 1 ^ 12 
or so-called Zeldovich-Frenkel equatio n 13 ' 14 for the one- 
component system. 

We start from the general Fokker-Planck equation for 
the multi-component system 1 ^ of the form 
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df(n) 
dt 
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FIG. 1. An example of the free-energy landscape of compos- 
ite nucleus that consists of a stable solid (n.2) surrounded by 
an intermediate metastable liquid (ni). This is calculated us- 
ing the capillarity approximation^ in the (metastable liquid 
(ni), stable solid (ri2))-space. The two red points indicate 
two saddle points with the same free-energy barriers which 
correspond to the critical nucleus of metastable liquid and 
the critical composite nucleus composed from solid core and 
surrounding metastable liquid wetting layer. The red line in- 
dicates the series nucleation route suggested by previous au- 
thors 2 ' 6 ' 7 . This free -energy landscape with two saddle points 
clearly indicates a parallel nucleation (two nucleation routes). 



where 



/(n) = /(ni,n 2 ,...) 



(2) 



is the number of clusters with composition n = 
(ril, n 2 , . . . ), and J is the nucleation flux vector 



J = J\&\ + J 2 e 2 



(3) 



where are unit vectors in the Cartesian (m, n 2 , . . . ) co- 
ordinate. By assuming the detailed balance at the equi- 
librium^, the flux components J, are derived from the 
potential $ defined by 



$(m,n 2 ,...) 



f(ni,n 2 , 



N (m, n 2 , . . . ) ' 



(4) 



where 



N{n 1 ,n 2 ,...) = JV expl — I (5) 

is the equilibrium number of clusters and W is the re- 
versible work (free-energy) of cluster formation with com- 
position m, TI2, ■ ■ ■ ■ 

Using the potential we can write the nucleation flux 

as 

Ji = -N (m, n 2 , • ■ • ) Ri (ni, n 2 , . . . ) — (6) 

where Ri represents the reaction rates of nucleation 
which will take different forms in vapor condensation^ 



and nucleation in condensed phased. Equations ([T]) 
and ([5]) lead to the Fokker-Planck equation^ or Frenkel- 
Zeldovic h 13 ' 14 equation when we consider only one com- 
ponent. Therefore, we have to investigate the Brownian 
motion on a multi-dimensional surface of the potential 

Since we are interested in the series and parallel nu- 
cleation phenomena, which are the most basic elements 
of complex nucleation, we consider the nucleation by a 
two-dimensional potential $ (711 , n 2 ), or the binary nucle- 
ation of two-component system that has been considered 
thoroughly^— . We concentrate on the steady-state nu- 
cleation rate only though previous authors^— were in- 
terested mostly in the relations between the nucleation 
current, the steepest-descent direction of the free energy 
landscape (work of formation) and the gradient direction 
of the potential $. 



A. Parallel nucleation 

In the two-dimensional space (ni, n 2 ), we have to solve 
the Fokker-Planck equation: 



df(n) 



V • J 



dJi dJ 2 



dt dm dn,2 

subject to the boundary condition^ 

/(|n|->oo) = 0, f (\n\ ^ 0) = N (n) 



(7) 



(8) 



The steady-state nucleation flux J st can be studied when 
df /dt = 0, which leads to 



<9Jf <9JJ 



dni dn-2 



= 



(9) 



from Eq. ([?])• Since divJ st = 0, in arbitrary closed re- 
gion on the two-dimensional space (m, n 2 ), the flux must 
be conserved. In other words, the flux coming into this 
region must be balanced to the flux going out from this 
region. 

It is customary to assume the nucleation flux flow 
through the saddle point of the free-energy landscape 
W (ni, n 2 )^I~— . Suppose we have two saddle points, 
say a and b whose positions are almost in parallel along 
the line ri\ + n 2 ~ constant. Therefore, the two saddle 
points correspond to the two types of nucleus of almost 
the same size and with different composition as shown in 
Fig. [2ja). In such a case, it is reasonable to assume the 
nucleation flux J st will fork into two branches J^f that 
passes through the saddle point a and that passes 
trough b as shown schematically in Fig. [U[a), and the 
flux far from the two saddle points becomes the sum of 
the two fluxes 



JSt TSt I TSt 



(10) 



where these two fluxes satisfy the conservation law 

divJf=0, K = a,6 (11) 
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(b)A scenario of series nucleation 



FIG. 2. A typical scenario of (a) parallel and (b) series 
nucleation on a free-energy landscape of the simplest two- 
component system: The nucleation flux J st indicated by the 
wide arrow must be continuous and conserved as there is no 
source term in Eqs. ([9]) and In (b), the problem is re- 

duced to the one-dimensional problem along the fictitious co- 
ordinate m. 



However, the steepest-descent direction of neither the 
free-energy landscape^ W nor the potential landscape 
$ does not necessarily tell you the direction of nucle- 
ation flux since the reaction rate Ri can be different 17 

In order to take into account the anisotropy of R, ac- 
cording to Trinkaus 18 and Wilcmsk: 22 , we introduce new 
variables v K = (v K ^,v K p) through 

^-< 1 = « 1 ) 1/ % M , (15) 

n 2 -n* K j = (R* K>2 ) 1/2 v K , 2 , (16) 

at each saddle point n = a,b, where i?* i are the reaction 
rates at the two saddle points. Then the free energy 
Eq. (fl~2|) takes the form 

W[nun 2 )-W; = \Y, y ^ti3 v ^ K = a > b 
where 

K,a = (Ki) 1/2 w* K>ij (R* Kd ) 1/2 (18) 

includes the information of not only the free energy land- 
scape W but also the reaction rates R. By diagonalizing 
Eq. HBD as 

W-W* = \(\ K e K + l K r? K ) (19) 

using the orthogonal transformation 22 

= v K ,l cosa K + v Ka su\a K (20) 
T) K = -^ K ,isina K + ^ K , 2 cosa K (21) 

with 



independently around the two saddle points a and b. This 
simple equation (fT0|) which represents the superposition 
of two independent nucleations has already been sug- 
gested and used by Wagner and StreyiS to explain their 
experimental result of nucleation in water/ n-nonane va- 
por mixture. 

General solutions for the steady-state nucleation rate 
at the saddle point were given by Trinkaus^. First of 
all, we expand the work of formation W (n) around the 
saddle points n* = (n* 1; n* 2 ) , (k = a, b): 

W (m, na) a W* + - J2 (ni ~ < % ) W* tij (n, - n%J) , 

K = a,b (12) 

where 

W: = W(n* KA ,n* K<2 ) (13) 
is the free-energy barrier at the saddle point k = a,b and 
/ d 2 W \ 

W*^ = U Z ■ ( 14 ) 

\ ' 3 / n— rt* 



tana K = (r„, u - r K , 22 - G K ) / (2T Ktl2 ) (22) 

and 

g k = V / ( r ^ii-r^ 22 ) 2 + 4(r K , 12 ) 2 , (23) 

we can obtain the eigenvalues X K and 7 K given by 

A K = (r K41 +r^ 22 -G K )/2<0, (24) 
7« = (r«,n+r re , 22 + G K )/2>0. (25) 

Then, the magnitude of the steady-state nucleation rate 
J® at the saddle point k = a, b is given simply by 



using the eigenvalues A K and 7 K 18 ' 22 . The total nucle- 
ation rate can be calculated by Eqs. (ITU1) and (|26l) . There- 
fore, the total nucleation rate will be larger than the in- 
dividual nucleation rates J„* and as the extra channel 
of nucleation reaction exists. 
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J st and we have 



FIG. 3. Free-energy landscape W when there exist two saddle 
points m a and nit along the fictitious one-dimensional axis m 
shown schematically in Fig. [2jb) . 



Since N* = N exp (-W*/kT), the two fluxes Jf and 
Jp will contribute almost equally when the free-energy 
barriers W* are almost the same magnitude. An example 
of the free-energy landscape when two saddle points have 
the same free-energy barriers has already be shown in 
Fig.ffl 

Even when two competitive routes k = a and n = b do 
not have the same free energy barriers, the corresponding 
two competitive nuclei nucleate probabilistically. A qual- 
itatively similar conclusion was reached by Sanders et 
alj^i using the artificially designed ^states Potts model. 
The probability is in fact proportional to the Boltzmann 
factor N* oc exp {-W*/kT) from Eq. 



B. Series nucleation 

When two saddle points a and b appear sequentially 
along the line which is almost perpendicular to the line 
ni + ri2 ~ constant, these two saddle points corresponds 
to the two types of nucleus of almost the same compo- 
sition n,2 oc n\ with different sizes (small and large). In 
such a case, it is reasonable to assume that the nucleation 
flux J will remain one stream as shown schematically in 
Fig. [2fb) and the problem will be reduced to the one- 
dimensional double-barrier— problem shown in Fig. [3] 

By changing the coordinate system along the flux, 
the problem can be mapped to the one dimensional 
problem with two barriers first studied by Valencia and 
Lipowsky 5 . Suppose we map the nucleation flux on the 
two-dimensional coordinates system {nx^n-i) to one di- 
mensional system with coordinate to (Fig. HJb)), then 
the nucleation flux Eq. ([5]) becomes one-dimensional 



dm 



f(m) 

R c ff (to) N (to) J dm \ N (to) 
f (to — > 0) / (m — > oo ) 



N (m 0) N (m ->• oo) 



(28) 



from Eq. (|27[) . Using the boundary condition^ in Eq. © , 
we found 



r 



dm 



o R c ff (to) N (m) 



(29) 



By dividing the integrals into two parts around the two 
saddle point m a and mb as 



J st = 



dm 



(m a ) Rett (m) N (TO) J( mb) R c ff (to) N (to) 



dm 



we arrive at the formula 

1 1 



(30) 



(31) 



for the series nucleation that is similar to the formula 
for the conductivity of two resistors connected in series, 
where 

J* =R eS {m K )N{m K )Z K , K = a,b (32) 



with the Zeldovich factor 



Z K = y/\d 2 W (m K ) /dm 2 \ jlifhT 



(33) 



is the nucleation rate when the saddle point n = a,b 
exists independently. The factor l/y/2-KkT appears 
in Eq. (|33[) of one-dimensional problem in contrast to 
Eq. (|26[) of two-dimensional problem since the general 
formula for the nucleation rate by Trinkaus^ for the d- 
dimensional free-energy surface is proportional to a factor 
(l/27rfcT) 1_d/2 . Eq. §T§ was initially derived by Valen- 
cia and Lipowsky^ using the Master equation and later 
re-derived by Valencia^ using the Kramers theory. Here, 
we have directly used the continuum description and the 
Fokker-Planck equation. Reduction of the nucleation 
rate J s from an independent J* is due to the escape flow 
from intermediate well (Fig. [3]! 26 . The double-barrier 
problem was also studied using Kramers theory by Nico- 
lis and Nicolis^I. Since they did not use the standard 
theory of nucleation based on the Fokker-Planck or the 
Zeldovich-Frenkel equation, they could not derive a sim- 
ple expression for the total nucleation rate like Eq. ([3H . 



("»s: (s^)- < 27) 

where we have used the definition of the potential $ in 
Eq. (HJ) and introduced an effective reaction rate R e ff (to). 
In the steady state, the nucleation flux is constant J = 



III. CONCLUDING REMARK 

In this report, we have summarized the basic ingre- 
dients of parallel and series nucleation when there exist 
two saddle points on the free-energy landscape. Exten- 
sion to the more complex nucleation process is possible 
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FIG. 4. Complex nucleation process represented by an anal- 
ogy to the electric circuit when there exist multiple saddle 
points a, b and c on the free-energy landscape W . Total 
steady-state nucleation rate J st is given by 



i 



by drawing the diagram similar to the electric circuitry 
as shown in Fig. 0J Figure 4 shows an example of com- 
plex nucleation when there are three saddle points on 
the free energy landscape. In this case, the steady-state 
nucleation rate will be calculated from 



1 



1 



1 



Jt + J? J? ' 



(34) 



Extension to more complex situation is obvious. 

However, it must be noted, our analysis is completely 
confined to the steady-state process. Therefore, transient 
properties are out of our scope, which can only be stud- 
ied numerically by solving coupled Master equatio n 20 ' 28 . 
Also, our analysis assumed that the nucleation flux goes 
through the saddle point. Saddle point avoidanc o 18 ' 21 
will be important if the anisotropy of the reaction ma- 
trix R is large or the ridge between the saddle points is 
law, which can occur at high temperatures or near the 
spinodal point. In such a case, ridge-crossing rather than 
saddle-crossing may occur. Then the nucleation flux will 



spread over the whole phase space and our description 
may break down. 
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